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ABSTRACT 


This NASA Grant Final Report is divided into two (2) 
parts. The first part summarizes the main effort of the grant- 
to investigate and develop precession maneuver control laws 
for single-spin spacecraft so that nutation is concurrently con 
trolled. Analysis has led to the development of two (2) types 
of control laws employing precession modulation for concurrent 
nutation control. Results have been verified through digital 
simulation of a Synchronous Meteorological Satellite (SMS) con- 
figuration . 

In the second part, an addition research effort, not 
originally considered, was undertaken to investigate the 
cause and elimination of nutation anomalies in dual-spin 
spacecraft. A literature search has been conducted and a 
dual-spin configuration has been simulated to verify that 
nutational anomalies are not predicted by the existing non- 
linear model. At the termination of the grant, the dual- 
spin research is still preliminary and no conclusions can be 
drawn as to the cause of the observed nutational anomalies in 
dual-spin spacecraft. 
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PART I. NUTATION CONTROL IN SINGLE-SPIN SPACECRAFT 


The purpose of Grant NGR 26-003-069 was to investigate 
variations in the basic precession maneuver control laws 
for single-spin spacecraft so that nutation is concurrently 
controlled. Nutation induced by precession was analyzed in 
detail. Subsequently, two control laws employing precession 
modulation for concurrent nutation control were developed 
and verified through digital simulation of a Synchronous 
Meteorological Satellite (SMS) configuration. Results were 
reported in the Semi-Annual Status Report "Nutation Control 
During Precession of a Spin-Stabilized Spacecraft" in December 
1973. 

Manuscripts describing the research performed in this part 
of the grant effort have been submitted to the American In - 
stitute of Aeronautics and Astronautics Journal of Spacecraft 
and Rockets in their Synoptic/backup paper format. Appendix 
A contains the Synoptic manuscript and Appendix B contains 
the backup paper manuscript, both entitled, "The Use of Pre- 
cession Modulation for Nutation Control in Spin-Stabilized 
Spacecraft" . 

The grant to the University of Missouri - Rolla, through 
direct financial aid or research involvement, supported the 
obtainment of the Master of Science in Computer Science degree 
by Mr. Vehbi Tasar and the Master of Science in Electrical 
Engineering degree by Mr. Richard J. Donner . 
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PART II. NUTATION ANOMALIES IN DUAL-SPIN SPACECRAFT 

As an additional research effort, not originally con- 
sidered in Grant NGR 26-003-069, a preliminary study was 
undertaken to investigate the cause and elimination of 
nutation anomalies in dual-spin spacecraft. The intent was 
to look at anomalous nutation behavior in dual-spin space- 
craft, in general, and, in particular, with regard to three 
NASA spacecraft, ITOS-D, ITOS-F, and AE-C. This behavior 
is characterized by either the decay or growth of induced 
nutation to some non-zero terminal value. 

A literature search was undertaken and a partial bib- 
liography is included in Appendix C. The published litera- 
ture on dual-spin dynamics generally fails to provide in- 
sight into mechanisms which might be responsible for this 
anomalous nutation behavior. However, [16] describes an 
effort to determine the cause of apparently similar behavior 
for a particular dual-spin spacecraft. Mechanical cross- 
coupling due to dissipation in the bearing assembly connecting 
the two spinning bodies was identified in this case as the 
probable cause. 

Yet another study [20] pin pointed the cross-coupling 
mechanism as electronic in nature. Here nutation influenced 
the position sensor and produced errors in the despin control 
torque. 

Communications with NASA regarding ITOS-D, ITOS-F, and 
AE-C have also suggested the interface between the spinning 
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bodies as the area wherein anomalous nutation behavior is 
generated. Prior to investigating the interface model, the 
dual-spin configuration described by Phillips [20] was 
simulated to verify that nutational anomalies were not pre- 
dicted by the existing nonlinear model. No nutational 
anomalies were observed in the simulation of the dual-spin 
model of [20], 

At the termination of Grant NGR 26-003-069, the dual- 
spin research was still in a preliminary stage and no conclu- 
sions could be drawn as to the cause of the observed nutational 
anomalies in dual-spin spacecraft. 
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APPENDIX A 

AIAA-JSR SYNOPTIC — THE USE OF PRECESSION MODULATION FOR 
NUTATION CONTROL IN SPIN-STABILIZED SPACECRAFT 
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NOMENCLATURE 

- accelerometer gain 

= accelerometer output 

= moments of inertia about principal axes 
= torque 

= spacecraft spin period 

- spacecraft nutation period 

= -1/x + ift, complex energy dissipation 

- angular position of control thruster 
with respect to spacecraft axes. 

= thrust pulse duration, in angle of 
nutation cycle 
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= thrust pulse duration, in seconds 

= accelerometer angular position 

= nutation phase angle at initiation of 
thrust pulse 

= measure of moment of inertia ratio 
between transverse and spin axes 

= energy dissipation time constant 

= nutation angle 

= transverse angular velocity 

= spacecraft spin rate 

= nutation frequency 

THEME 


Spin stabilization is used to inertially fix a spacecraft's 
thrust vector during a transfer ellipse; however, spin about an 
axis other than that of maximum inertia results in an unstable 
equilibrium condition. Energy dissipation resulting from fluid 
motion in heat pipes and fuel tanks causes the nutation angle to 
increase. Furthermore, during the transfer ellipse, if it is 
necessary to precess the spacecraft's spin axis to reorient the 
thrust vector, this maneuver can also cause nutation to increase. 

This paper analytically derives relations which determine the 
nutation effects induced in a spinning spacecraft by periodic 
precession thrust pulses. By utilizing the idea that nutation need 
only be observed just before each precession thrust pulse, a tedious 
continuous time derivation reduces to a simple discrete-time derivation 
which can be easily solved through the use of z-transf orms . 

The analytic results obtained are used to develop two types of 
modulated precession control laws which use the precession maneuver 
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to concurrently control nutation. 

CONTENTS 


In 1968 Grasshoff [1] published the original work on a control 
law which can be used to sense and remove nutation automatically. 
Grasshoff 's basic differential equation for the transverse angular 
velocity of the spinning spacecraft is 

bi ( t) + ifia) (t) = N^e'*'^ (1) 


and the solution for a control thrust duration of (t^-t ) seconds 


is given by 


. -in (t -t ) . N o e 

u(t i) = “o e 1 0 + Tn 


i A 


(t l" t 0 ) ] (2) 


Grasshoff defines the nutation angle as the peak value of 


0 = (t) / (u^-fi) . (3) 

Nutation can be sensed by an accelerometer with sensitive axis 
parallel to the spacecraft's spin axis and can be removed by 
properly timing the firing of a small thruster whose thrust axis 
is also parallel to the spin axis. The sensed acceleration is 
given by 

a = a^Q sin (fit + y) . (4) 

Energy dissipation causes nutation to continually increase. 
For small pertubations , the increase in nutation angle, 0, is 
generally represented [3] as 

e ( t) = 0 o e t/T , (5) 

in which t is the time constant of energy dissipation. 
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Due to the relationship in (3) between 0 and w(t), energy 
dissipation can be included in (1) to give 

Jj(t) + (-i + ifl)w(t) - No elX 

Precession of the spin axis of a spinning spacecraft can be 
effected by using the same thruster that is required for nutation 
control. For a precession maneuver consisting of a thrusting sequence 
of equally spaced thrust pulses that are At wide and have a period 
= 2 it/u> s , (1) becomes 

N, n-T <t<n*T + At, 

s- s 

n=0 ,1,2... ( 7 ) 

O , n-T +At< t< (n+1) T . 

s s 

The continuous solution for (7) is complicated by the 
essentially infinite pulse train. The inverse Laplace transform 
of the Laplace transform of (7) is difficult because the pole at 
s = -ift cannot be separated by a contour from the poles at ±i2mr/T s , 
n-1,2,3,.... (See reference (4] for a discussion of Laplace trans- 
forms of pulse trains.) 

The solution of (7) is straightforward if it is sufficient 
to know only the result at discrete times, for example, at the dis- 
crete point, k T , just prior to the next thrust pulse. Thus, (7) 
can be changed to a first-order difference equation. (See reference 
[5] for a discussion of difference equations and z-transform theory.) 



By taking the z-transform and then taking the inverse z-transform, 

the discrete time solution for (1) is 

r -ift (nT ), N 0 e ,, i«At) -ifiT , 1 , e~ in(nT s ) , 

U [n T]=u e s' +nr (l-e e s[ _ lflT + — 5 <8) 

1-e s 1-e s 

The discrete time solution for (6) can be similarly obtained as 
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N rt e^ . i -a (nT ) 

r m , -a (nT ).0 .. -aAt)_-aT_ r 1 , e s 

W nT s ^o e S a e ^ -aT aT 

* u 1-e s 1-e s 


(9) 


where a = -(1 /t) + in. 

The conclusions are that with no energy dissipation and with 
a judicious choice of inertia ratios, which determine the relation- 
ship between T and T n , nutation resulting from unmodulated pre- 
cession is bounded and periodic and can be kept at tolerable levels. 
With energy dissipation, nutation resulting from unmodulated pre- 
cession is unbounded. Thus, some sort of nutation control is required 
during a precession maneuver. 

Because precession and nutation control can use the same 
thruster, precession thrust pulses can be selected that will con- 
currently reduce nutation. The effect of a single precession pulse 
on nutation can be demonstrated by a graphic representation of (2) 
shown in Figure la. 

The range of the initial position 0 Q can be analytically derived 
from (2). After trigonometric manipulation, replacement of w by 
-i4>o, and (t^-t Q ) by AO, the expression for the range of 0 Q , for 
which o)(t^)<w(t ) is given by 

arccos [-m sin (AO/2) ] ^-A - (AO/2 ) <0^< 

arccos[-m sin (AO/2 )] ^.^-A- (AO/2 ) , (10) 


in which arccos [ • ] XI and arccos ( - ] j XI refer to angles in the II 
and III quadrants, respectively, and m is given by 

m = (N /fl) (1 /oj (t ) ) . (11) 

To summarize, (10) shows that a relation exists which can be 
used to modulate, or gate, the precession pulse train so that only 
precession pulses are allowed which concurrently reduce nutation. 
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In fact, the nutation sensing accelerometer discussed previously 
and a threshold derived from (10) can be used to provide precession 
pulse gating. 

Figure lb shows the use of the accelerometer output for time- 
optimal nutation control. Any thrust which has a pulse duration of 
less than tt/^ but occurs during the negative half-cycle of the 
accelerometer output, will also reduce nutation, but suboptimally . 

Figure lc is a block diagram of a modulated precession, nuta- 
tion control. The end result is that only precession pulses that 
concurrently reduce nutation are passed to the thruster. 

Two precession modulation nutation control methods are now 
presented. The first provides continuous precession modulation, 
and the second allows unmodulated precession until nutation increases 
to a threshold value at which time the precession is modulated. The 
corresponding figures show digital computer simulation results of 
these control laws applied to the Synchronous Meterological Satellite 

(SMS) configuration specified as follows: mass properties: mass = 

2 

42.8 slugs, I - 244.4, I = 246.9, I - 97.3, all in slug-feet , 
x y 2 

c = 0.6039? angular velocity: w = 90 rpm, Q = 5.69 rad/sec; 

b 

thruster properties: thrust = 5 lbs., moment arm = 3 feet, duty 

cycle for precession = 1/12 second (45° of spin period) , duty cycle 
for nutation = 0.552 seconds (1/2 nutation period); energy dissipation 
time constant: t = 180 seconds. 

Figure 2a shows the behavior for the continuous precession modu- 
lation nutation control. Because the gate width, as given in (10) , 
varies inversely with nutation, the precession increase and nutation 
decrease will be fairly rapid initially. As the nutation reduces. 
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fewer precession pulses are passed, and the increase in precession 
with respect to time reduces to the constant rate shown. Consequently, 
as shown in Figure 2a, 2000 seconds are required to precess 30 , but 
after the first ten seconds of the precession maneuver nutation remains 
below 2 milliradians . 

Figure 2b shows the behavior for the nutation control law which 
allows unmodulated precession until nutation increases to a thres- 
hold value of 5.48 milliradians at which time the precession is modu- 
lated. Thus, the precession and nutation increase until the nutation 
threshold is reached. At this time, precession modulation occurs. 

The slope of the precession increase diminishes as the result of 
precession modulation, and the nutation is held near the threshold 
level . 

In the second case the control law passes all precession pulses 
until nutation builds up to a threshold of 5.48 milliradians. Conse- 
quently, precession is more rapid. As figure 2b shows, 2100 seconds 
are required to precess the entire 180°. 

By way of contrast, a simulation of a conventional control law 
which allows nutation control, or precession, but not concurrently, 
for the SMS configuration described requires 1990 seconds to precess 
180°. The fuel budget for the modulated precession control law shown 
in Figure 2b is 228.9 seconds of thruster fuel. For the conventional 
type of precession and nutation control the fuel budget is 307.9 seconds 
of thruster fuel. 

In conclusion, for spacecraft such as those described by Grass- 
hoff [1], and Taylor [2], for the SMS, and others emphasizing active 
on-board nutation control with an accelerometer as the nutation sensing 
element, the mechanization of the control laws which use the precession 
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maneuver to concurrently control nutation is straightforward and 
requires little additional electronic circuitry. 
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Figure 1. 


Figure 2. 


(a) Graphic Representation of Equation (2) . 

(b) Nutation Sensing Accelerometer Output. 

(c) Modulated Precession Nutation Control. 

(a) SMS Simulation For Full Modulated Pre- 
cession Nutation Control. 

(b) SMS Simulation For Full-Precession-to- 
Threshold Modulated Precession Nutation 
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NOMENCLATURE 


= accelerometer gain, Eq . (5) 

= accelerometer output, Eq - (5) 

ss moments of inertia about principal axes 

= least common multiple of T , T^ 

= torque, Eqs. (1), (3) 

= torque, Eqs. (2), (3) 

= spacecraft spin period 

= spacecraft nutation period 

= -1/t + ift, complex energy dissipation, 

Eq. ( 14 ) 

= angular position of control thruster with 
respect to spacecraft axes. 

= thrust pulse duration, in angle of 
nutation cycle 

= thrust pulse duration, in seconds 

= accelerometer angular position, Eq. (5) 

= nutation phase angle at initiation of 
thrust pulse 

= measure of moment of inertia ratio 
between transverse and spin axes 

= energy dissipation time constant 

— nutation angle, Eqs. (4), (6) 

= transverse angular velocity 

= spacecraft spin rate 

= nutation frequency 
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I . INTRODUCTION 

Spin stabilization is used to inertially fix a spacecraft’s 
thrust vector during a transfer ellipse; however, spin about an 
axis other than that of maximum inertia results in an unstable 
equilibrium condition. Energy dissipation resulting from fluid 
motion in heat pipes and fuel tanks causes the nutation angle to 
increase. Furthermore, during the transfer ellipse, if it is 
necessary to precess the spacecraft’s spin axis to reorient the 
thrust vector, this maneuver can also cause the nutation to increase. 

This paper analytically derives the relations which determine 
the nutation effects induced in a spinning spacecraft by periodic 
precession thrust pulses. By utilizing the idea that nutation need 
only be observed just before each precession thrust pulse, a 
tedious, unmanageable continuous time derivation reduces to a simple 
discrete-time derivation which can be easily solved through the use 
of z-transf orms . 

The analytic results obtained are used to develop two types of 
modulated precession control laws which use the precession maneuver 
to concurrently control nutation. 

Section II briefly reviews nutation in spinning spacecraft and 
active nutation control. Section III develops the relations which 
determine the nutation effects induced by periodic precession thrust 
pulses. Section IV develops the relations required to modulate the 
precession thrust pulse sequence so that nutation is reduced. Sec- 
tion V describes two types of modulated precession nutation control 
laws which are demonstrated in the digital simulation of a Synchronous 
Meteorological Satellite (SMS) configuration. Finally, Section' VI 
presents conclusions. 
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II. NUTATION AND ACTIVE NUTATION CONTROL 

In 1968 Grasshoff [1] published the original work on a 
control law which can be used to sense and remove nutation auto- 
matically. Grasshoff* s basic differential equation for the trans- 
verse angular velocity of the spinning spacecraft is 

w(t)+ifta>(t)=N (1) 


and the solution for a control thrust duration of (t^-t ) 
is given by 


N e 

x -ia(t,-t_) ^ o 

“'V = U) 0 e 10+^ 


iX 


[l-e _ifl(t rV] 


seconds 


( 2 ) 


where N in (1) is defined as 

N = N q e iX (3) 

and X is the angular position of the thruster in the spacecraft’s 
coordinate system. Thus, if a thrust pulse is initiated at the 
proper time and for the proper duration, w(t^) is less than a)(tg) f 
and the transverse angular velocity is reduced. 

Grasshoff proceeds to define the nutation angle as 


0 - w(t)/(io -fl). (4) 

Nutation is sensed by an accelerometer with sensitive axis 
parallel to the spacecraft’s spin axis and is removed by properly 
timing the firing of a small thruster whose thrust axis is also 
parallel to the spin axis. 

The sensed acceleration is given by 
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a = a Q 6 sin (fit + y) (5) 

in which 6 is the nutation angle, y is the phase angle that 
defines the position of the accelerometer with respect to the 
position of the thruster, and a^ is a gain that is a function 
of the accelerometer’s radial position and the spacecraft's spin 
rate . 

Energy dissipation due to liquid movement in fuel tanks and 
heat pipes of the spinning spacecraft causes nutation to contin- 
ually increase. For small pertubations , the increase in nutation 
angle, 0, caused by energy dissipation is generally represented 
[3] as 

e(t) = 0 o e t/T , (6) 

in which x is the time constant of energy dissipation. 

Due to the relationship in (4) between G and w(t), energy 
dissipation can be included in (1) to give 

w(t) + (-- + ifi)w(t) = N (7) 

Figure 1 illustrates a simple, nonoptimal , nutation control 
law investigated by Taylor [2] for the International Ultraviolet 
Explorer (IEU) spacecraft. A threshold is established which rep- 
resents the amount of nutation that can be removed by a properly 
phased thrust pulse of duration equal to one spin period, i.e., 

2 Ti /oo . Nutation is sensed by the accelerometer. When the next 
s 

zero-crossing of the acceleration sinusoid occurs, the thruster 
fires and reduces the nutation to near zero. Illustrations la and 
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lb show the nutation control for a small initial nutation, and 
lc shows the series of repeated thrust pulses that are required 
for a large initial nutation. 

The nutation control law presently planned for the IUE 
operates in the same manner except that the thrust pulse is of 
a duration equal to one-half the nutation period, i.e., w/il. 

III. NUTATION DURING UNMODULATED PRECESSION 

Precession of the spin axis of a spinning spacecraft can be 
effected by using the same thruster that is required for nutation 
control; however, in this case, precession in a uniform direction 
is the result of a train of thrust pulses which have a frequency 
equal to the spin frequency of the spacecraft. This section shows 
how precession affects the nutation of a spacecraft with and with- 
out energy dissipation. 

Precession Without Energy Dissipation 

Equation (1) describes the change in transverse angular 
velocity, and by virtue of (6), the change in nutation, without 
energy dissipation, as a result of the forcing function N. 

Consider a precession maneuver consisting of a thrusting sequence 
of equally spaced thrust pulses that are At wide and have a 
period T = 2 tt/uj . Thus, 

N, n*T < t <n • T + At, 
s- s 

n=0 , 1 , 2 . . . 

0 , n*T +At< t< <n+l) T . 
s s 


a) (t) + iftu) { t) 


( 8 ) 
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The continuous solution for Eq. (8) is complicated, by the 
essentially infinite pulse train. For example, the Laplace 
transform of Eq. (8) is given by 


0). 


w ( s ) = 


s+iU 


-At s v 
N(l-e ) 

s(s+ift) (1-e ^s s ) 


(9) 


(See reference [4] for a discussion of the Laplace transforms of 
pulse trains.) The transient portion of Eq. (9) has poles at 
s=0, s = -ift, and s = +i2mr/T , n = 1, 2, 3.... 

b 

The inverse Laplace transform of Eq. (9) is extremely complex 
and tedious because the pole at s = -i^ lies on the imaginary axis 
and cannot be separated by a contour from the poles at + i2n7r/T s . 
Thus , a continuous time-domain solution is formidable and not worth 
the effort because a bound on Eq. (8) will suffice. 

The solution of Eq. (8) is straightforward if it is sufficient 
to know only the result at discrete times, for example, at the dis- 
crete point, k T , just prior to the next thrust pulse. Thus, Eq. 
(8) can be changed to the first-order, difference equation 


u>[k T g ] 


u>[ (k-l)T g ] 


-iQT 
e s 



. . Lil At 

(1-e 


% -iftT 
) e 


s . 


( 10 ) 


(See reference (5] for a discussion of difference equations and 
z-transform theory.) 

By taking the z-transform of Eq. (10) and then taking the inverse 
z-transform, the discrete time solution is 


r mi -i£2(nT ) 
a) [n T] =WQe s + 


N ,, iftAt) -iflT , 1 

15 (1 ' e e s[ : -ifl'f— 

1-e s 


-ift(nT ) 
e s 

, iftT 
1-e s 


3 (ID 
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Equation (11) can be simplified by combining constant para 
meters to 


ti) 


nT ] = 
s 


e 


-10 i 


r , » —i 0 i •, _ 10 2 

+ { [(0 +k n e ]e 
o 1 


}e 


-ifi(nT) 


( 12 ) 


in which 

N Sin (l/2Mt) ( 13 ) 

K 1 Sin (1/2S2T ) ' 

s 

and and 02 ate constant angles. Equation (12) is periodic 
with period, £, where £ is the least common multiple of [T g , T fi ] , 
and T^ is the nutation period 2tt/Q. In Eq . (13) , is bounded 

except for the singularity that occurs when T g =Tr/ft, which is the 
case when the spin period equals one-half the nutation period. 
Thus, nutation resulting from precession, in the absence of energy 
dissipation, is bounded except at one singularity. 


Precession With Energy Dissipation 


Equation (7) describes the change in transverse angular 
velocity, and consequently nutation, with energy dissipation. 
As in the case without energy dissipation, the use of the z- 
transform to provide a discrete-time solution offers the most 
straightforward approach. Thus, the discrete-time solution of 
Eq. (7) is 


a 1 


[nT ] s u e a(nT s ) + (N/a) (1-e aT s [ 


. -aT c 
1-e s 


-a (nT„ ) 
e 

1-e s 


(14) 


where a = -(l/i) + ift 
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Equation (14) can be simplified by combining constant parameters to 

( 15 ) 


r <nT s> -isMi.T > 

o) [nT ] - (u) -k ) e e “ s f k . 

sol 2 


in which w[nT ] is unbounded unless w exactly equals k , which 
s o 1 

is only theoretically possible. 

The conclusions are that with no energy dissipation and with 
a judicious choice of inertial ratios, which determine the rela- 
tionship between T g and , nutation resulting from unmodulated 
precession is bounded and periodic and can be kept at tolerable 
levels. With energy dissipation, nutation resulting from unmodu- 
lated precession is unbounded. 


IV. NUTATION DURING MODULATED PRECESSION 

Because precession and nutation control can use the same 
thruster, precession thrust pulses can be selected that will con- 
currently reduce nutation. 

The effect of a single precession pulse on nutation can be 
demonstrated by a graphic representation of Eq . (1) shown in 

Figure 2a. In this figure, AO is the duration of the thrust pulse 
with respect to the nutation cycle. Thus, for an initial w(t ) 
at an angle 0 q , under the influence of a thrust pulse of (t^-t ) 
seconds in duration, w(t^) is the resultant of the fixed thrust 
vector plus m(t ) rotated through the angle AO. For the conditions 
shown in Figure 2a, u>(t^) is greater than oi(t ), i.e., nutation has 
increased. 

Figure 2b shows the locus of u(t^) for a thrust pulse of given 
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duration as a function of o>(t ) and the initial angle 1 . Note 

o o 

that nutation is reduced whenever co(t^) occurs in the hatched 
area . 

Figure 2c is essentially the same as Figure 2b except that 

circle A and circle B are scaled to show the initial angle 

The displacement of the scales on circle A and circle B results 

from the value of A4>. Thus, for a thrust pulse of duration A4>=90°, 

starting when the tip of the u>(t Q ) vector is at <P q = 0° on circle 

B, the result will be an ^(t^) vector with the tip at 0° on circle 

A. Figure 2c also shows the range of for which the thrust pulse 

will reduce nutation, i.e., w(t, )<w(t ). In this case, 

l-o 


co (t n ) <co (t ) for 
1 - o 


0°<$ <29 

- o- 


239 ° < $ <360°. 

- o- 


The range of can be analytically derived from Eq. (2) . 

— i <J>o 

After trigonometric manipulation, replacement of , and 

(t,-t ) by A3> gives 
1 o 

-a-*, 2N - A . ■ , A<K 

, . , -l$ ft “lA$ , 0 . ,A1\ i(X - — sp) 

co(t 1 ) = L 0 Q e Oe + — sm(— )e 2 . (1 


The expression for the range of $, for which u> ( t^ ) <co ( t ) , is 
given by 

arccos [ -m sin(A$/2)J,. I -’A--(A$>/2)<$Q 

<arccos[-m sin ( A$/2 ) 3 ttt”X" ( A3>/2 ) , (18) 

in which arccos [•]■£-£ an d arccos [•] refer to angles in 

the II and III quadrants, respectively, and m is given by 


- (N o /ft) (l/oi <t Q ) ) . 


m 


(19) 
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Physically, m is the ratio of the radius of the thrust pulse 
lobe and the initial nutation. In Figure 2c, for example, 
m = 0.375. 

For a given precession maneuver, all variables in Eqs . (18) 

and (19) are fixed except for oMt^) which is proportional to the 
nutation, 8, at the beginning of a precession thrust pulse. 

To summarize, Eq . (18) shows that a relation exists which 

can be used to modulate, or gate, the precession pulse train so 
that only precession pulses are allowed which concurrently reduce 
nutation. In fact, the nutation sensing accelerometer discussed 
in Section II and a threshold derived from Eq . (18) can be used 

to provide precession pulse gating. 

Figure 3a shows the use of the accelerometer output for time- 
optimal nutation control. The duration of the nutation removal 
thrust pulse is it / ft, which is one-half the nutation period. Fur- 
thermore, any thrust, which has a pulse duration of less than ir /ft 
but occurs during the negative half-cycle of the accelerometer 
output, will also reduce nutation, but suboptimally . 

Figure 3b is a block diagram of a modulated, precession, 
nutation control. The end result is that only precession pulses 
that concurrently reduce nutation are passed to the thruster. 

To substantiate that a nutation control law of the type 
shown in Figure 3 will control nutation over an extended period 
of time, it is necessary to determine the average and maximum 
time between modulated precession pulses. Because precession 
pulses are used for nutation control, the nutation removed by a 
precession pulse must be greater than nutation buildup between 



28 


precession pulses resulting from energy dissipation. As nutation 

increases, the gate width in Figure 3b increases to 180°. Thus, 

on the average 50 percent of the precession pulses will be passed 

and 50 percent will be blocked. Nutation can be controlled if the 

nutation increase in time 2T , where T is the spin period, is 

s s 

less than the nutation removed by one precession thrust pulse. 

The maximum time between gated precession thrust pulses is 
difficult to show analytically and depends on the ratio of to 
T s which may be irrational. Computer solution for the spacecraft 
described in Table 1 shows that for a gate width of 50 percent, 
maximum precession pulse separation is 3T^ . Furthermore, for the 
spacecraft of Table 1, modulated precession pulses can maintain 
the nutation angle at less than 6 milliradians for precession pulse 
separation up to 50 to 60T s - 

V. MODULATED PRECESSION NUTATION CONTROL 

This section describes two control laws which use modulated 
precession pulses to control nutation. The corresponding figures 
show digital computer simulation results of these control laws 
applied to the Synchronous Meterological Satellite (SMS) configura- 
tion specified in Table 1. 

Of the two precession modulation nutation control methods 
described below, the first provides continuous precession modula- 
tion, and the second allows unmodulated precession until nutation 
increases to a threshold value at which time the precession is 
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modulated. Figure 4 shows the behavior for the continuous pre- 
cession modulation nutation control. Because the gate width, as 
given in (18) , varies inversely with nutation, the precession 
increase and nutation decrease will be fairly rapid initially. As 
the nutation reduces, fewer precession pulses are passed, and the 
increase in precession with respect to time reduces to the constant 
rate shown. 

Figure 5 shows the behavior for the nutation control law which 
allows unmodulated precession until nutation increases to a thres- 
hold value at which time the precession is modulated. Thus, the 
precession and nutation increase until the nutation threshold is 
reached. At this time, precession modulation occurs. The slope 
of the precession increase diminishes as the result of precession 
modulation, and the nutation is held near the threshold level. 


VI. CONCLUSION 

This paper analytically derives the relations which determine 
the nutation effects induced in a spinning spacecraft by periodic 
precession thrust pulses. By utilizing the idea that nutation 
need only be observed just before each precession thrust pulse, 
a tedious, unmanageable continuous time derivation reduces to a 
simple discrete-time derivation which can be easily solved through 
the use of z-trans forms . 

The analytic results obtained are used to develop two types 
of modulated precession control laws which use the precession 
maneuver to concurrently control nutation. In the first case the 
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control law passes only those precession pulses that concurrently 
reduce nutation. Consequently, as shown in Figure 4, 2000 seconds 
are required to preccss 30°, but after the first ton seconds of 
the precession maneuver nutation remains below 2 mi 11. i radians . 

In the second case the control law passes all procession 
pulses until nutation builds up to a threshold of 5.48 milli- 
radians. For nutation above the threshold only nutation reducing 
precession pulses are passed. 

By way of contrast, a simulation of a conventional control 
law which allows nutation control or precession, but not con- 
currently, for the SMS configuration described in Table 1 requires 
1990 seconds to precess 180°. The fuel budget for the modulated 
precession control law shown in Figure 5 is 228.9 seconds of 
thruster fuel. For the conventional type of precession and nuta- 
tion control the fuel budget is 307.9 seconds of thruster fuel. 

In conclusion, for spacecraft such as those described by 
Grasshoff [1], and Taylor [2], for the SMS, and others emphasizing 
active on-board nutation control with an accelerometer as the 
nutation sensing element, the mechanization of the control laws 
described herein is straightforward and requires little additional 
electronic circuitry. 
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TABLE I 


SMS SPACECRAFT CONFIGURATION 


Mass Properties 

Mass - 42.8 Slugs 

I - 244.4 Slug-Feet 2 

x 

2 

I = 246.9 Slug-Feet 

y 

I = 97.3 Slug-Feet 2 

a = 0.6039 

Angular Velocity 

to = 90 RPM (spin axis angular velocity) 

ft = 5.6916 rad/sec {nutation frequency) 

Thruster Properties 

Thrust = 5 lbs. 

Moment Arm = 3 Feet 

Duty Cycle = 1/12 seconds (45° of spin period) for Precession 
= 0.552 seconds (1/2 nutation period) for nutation 

Energy Dissipation Time Constant 
t = 180 seconds. 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


FIGURE CAPTIONS 


(a) Nutation Control for Small Angles (Phase Plano) . 

(b) Nutation Control for Small Angles. 

(c) Nutation Control for Large Angles. 

(a) Graphical Interpretation of Equation (1) . 

(b) Nutation Reduction Area 

(c) Nutation Reduction Area Scaled to Initial Angle 

(a) Nutation Sensing Accelerometer Output 

(b) Modulated Precession-Nutation Control 

SMS Simulation Results for Full Modulated Precession 
Nutation Control 


SMS Simulation Results for Full Precession/Threshold 
Modulated Precession Nutation Control 
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